Plastic response of dislocation glide in solid helium under dc strain rate loading 
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In this work we develop a model for the gliding of dislocations and plasticity in solid ''He. This 
model takes into account dislocation multiplication, Peierls barrier, and dislocation interactions. We 
specifically examine stress-strain curves and temperature effects. As a function of temperature, we 
observe plastic deformation and discuss how this may be related to the experimental observation 
of elastic anomalies in solid hep ''He that have been discussed in connection with the possibility of 
supersolidity. Our model predictions can be compared directly to constant strain rate experiments 
and provide bounds on model parameters. 
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I. INTRODUCTION 

There has been recent excitement about the observa- 
tion of a softening effect of the shear modulus in solid 
hep ^He as a function of increased temperature, which 
was interpreted in terms of damped vibrating dislocations 
in a dislocation networkii"— Since the anomaly in *He 
correlated with the previously reported supersolidity, as 
seen by a period drop in torsional oscillator experiments 
in the same temperature range,—"— it was argued to be 
closely related to the same phenomenon and provide ev- 
idence for dislocations playing an important role in the 
supersolid behavior of a Bose solidjii^i^ More recently, 
there have been more precise follow-up experiments on 
solid ^He in porous Vycor glass, where the period drop 
is no longer observed)^ or torsional oscillator measure- 
ments of shear-strained ^He, where the period drop and 
accompanying dissipation could be drastically enhanced 
at high temperatures, indicative of work softening and a 
nonsuperfluid originJ^ A still open question is what is 
responsible for the anomalies in the shear modulus, if it 
is not caused by an interplay between supersolidity and 
dislocation motion,^- 

One possible cause for the change of elastic proper- 
ties measured in experiments could be plastic deforma- 
tion of the type that arise from the gliding motion of 
dislocations, in contrast to the vibrating string model 
of dislocations pinned by a network, which so far was 
used primarily to interpret experiments.—"— If disloca- 
tion glide occurs in solid ^He, there should be clear sig- 
natures in experimentally measurable quantities such as 
stress-strain curves. This connection is well established 
for other materials, such as metals, and discussed in 
the literature.—"— The theory of dislocation motion and 
plastic deformation of crystals has a long tradition and 
one typically studies the temperature and strain rate ef- 
fects of the stress-strain curve of the crystal .^°'^^" — For 
that purpose a model based on rate-controlled mecha- 



nisms is necessary to describe the dynamics of disloca- 
tions, where thermally-activated flow is assumed with 
basal or prismatic slip in hep crystals (like solid ^He at 
low temperatures) . In addition, the theory must account 
for the creation and removal of dislocations, as well as 
their glide. Since mobile dislocations get pinned many 
times and at many points along their length due to the 
interactions with other dislocations, the effect of dislo- 
cation interaction must be included specifically. Finally, 
the rate of plastic strain is controlled by the rate at which 
thermal energy assists these segments to overcome their 
energy barriers, allowing the rest of the dislocation line to 
bow out and spread, before it is stopped agaiui ^^i^^i^^" — 

Many of the recent studies of the elastic properties 
of solid ^He were performed with ac measurement tech- 
niques of alternating shear strain, where the piezoelectric 
transducers were embedded in the solid helium. How- 
ever such ac measurements can present certain difficul- 
ties, because the boundary conditions, loading methods, 
and strain fields at the corners of the transducers can 
have a significant influence on the measured materials 
properties in solid *He<^ In contrast, here we propose 
a controlled dc drive experiment at constant strain rate 
with structural heterogeneities to examine the effects on 
the stress-strain curve. Our model results suggest that 
a dc drive is probably more appropriate for establish- 
ing whether dislocation glide is occurring versus the vi- 
brational motion of network pinned dislocations,—"—!^ 
since it avoids complications due to pile-up and annihila- 
tion of dislocations moving back and forth as can occur 
in an ac setup. 

In the theory of plastic deformation of crystals, the 
role of dislocations gliding in specific crystallographic di- 
rections on specific crystallographic slip planes is related 
to a change of the work-hardening rate of the solid with 
temperature and strain rate i^^'^^" — The hypothesis that 
this microscopic picture of slip and propagation of dis- 
locations is relevant to solid helium received additional 



support from recent experiments in the hep phase of solid 
"'Hej^ Naturally, the same mechanism may also be re- 
sponsible for similar anomalous behavior in disordered 
bcc crystals of "'He at temperatures close to the melting 
line^dL and in hep ^He,^ 

In this paper, we develop a dislocation glide model for 
solid "'He that is motivated by earlier work of Kocks and 
coworkers2£i2^iH on the dynamics of mobile dislocations 
moving through a dislocation forest. This model is based 
on the thermodynamics and kinetics of dislocation glide 
to describe plastic deformation of crystals. However, we 
extend the framework of this classical plasticity model 
by interpreting the zero-point motion and quantum ef- 
fects of solid helium as not qualitatively changing the 
kinetics and dynamics, but instead effectively renormal- 
izing certain thermodynamic quantities, such as by low- 
ering the pinning potential, Peierls barrier, or screening 
of the dislocation-dislocation interaction42r— In the case 
of dc shear strain rate experiment, e = const., with shear 
stress T, we predict the softening of the work-hardening 
rate (WHR). dr/de, with increasing temperature T due 
to dislocation glide and the creation and multiplication 
of dislocation loops. In the limit of e ^ the WHR 
approaches the elastic shear modulus, dr/de — > ^, and 
exhibits purely elastic behavior. Our dislocation glide 
model also predicts that higher "^He impurity concen- 
tration in solid ^He pushes the onset of the magnitude 
change in the WHR to higher temperatures, assuming 
that more ^He atoms equals more pinned dislocations. 
Pinning increases the number of immobilized dislocations 
and causes the crystal to harden, since no plastic defor- 
mation work is performed by them. At the same time, the 
zero-temperature value of the WHR is unaffected by im- 
mobile dislocations, whereas its high-temperature value 
is increased with their increased density — the crystal be- 
comes hardened through work performed on it. Finally, 
our calculations show that at zero temperature the WHR 
equals the elastic shear modulus independent of applied 
strain, while at finite temperature it decreases with in- 
creasing strain to a saturated value4^ 

Since we explicitly incorporate the effects of plastic 
deformation and interactions between dislocations, our 
dynamic glide model goes beyond the theory of vibrat- 
ing dislocations, which has been successful in describing 
the damping of pinned dislocations with small amplitude 
vibrationsji^rJ^iM Furthermore, in the spirit of recent 
approaches including quantunt^ and thermal^ fluctu- 
ations for the process of overcoming the Peierls barrier, 
we account for fluctuations in the thermodynamic part 
of the model with renormalized phcnomenological model 
parameters of the energy barrier. A key result of this 
work is the prediction of a dissipation peak caused by 
plastic deformation of the solid in the same tempera- 
ture region, where the WHR changes most rapidly. The 
similarity of the dc strain rate model predictions with 
existing shear measurements of applied ac strain rate 
is remarkable^^i^^— although a direct comparison is not 
possible due to significantly different loading and revcrse- 




FIG. 1: (color online). Characteristic dislocation network in 
crystalline materials illustrates the complexity of dislocation 
glide, (a) Snapshot of simulation of typical three-dimensional 
(3D) network of dislocations in crystal. Different colors repre- 
sent different Burgers vectors that indicate the magnitude and 
direction of strain imparted by glide of the dislocations, (b) 
Top view of 3D dislocation structure. Pinning points marked 
by arrows will unzip from junctions under external loading, 
see the SM for details on the simulation and the animation4i 



loading histories. Therefore, our results for the proposed 
dc strain rate experiments will provide stringent tests for 
the study of the plastic response of dislocation glide in 
the classical versus quantum regime at low temperatures 
of solid *He. 



II. PLASTIC GLIDE OF DISLOCATIONS 

In crystalline materials like "^Hc, plastic deformation 
can arise from glide of linear structural defects called dis- 
locations. Dislocations are commonly found in "'He crys- 
tal, even in the cleanest samples^ due to the growth and 
cooling processes. A typical dislocation structure simi- 
lar to that in metals is assumed in our model. In Fig. [1] 
wc show a snapshot of the complexity of such a three- 
dimensional (3D) dislocation network commonly studied 
in metals, see the Supplemental Material (SM)^ for de- 
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tails on the simulation and the animation of the dislo- 
cation glide, as well as the unzipping of pinning points 
marked by arrows. We borrow various elements from the 
seminal works on evolution laws of thermally activated 
dislocation glide by Essmann and Mughrabi)^ and Kocks 
and collaborators and categorize the dislocations as 
being either stored or mobile. The stored dislocations can 
assemble into networks inside individual crystals and do 
not contribute to the plastic deformation. Hence an in- 
crease in pinned or stored dislocations increases the WHR 
of the solid at higher temperatures. The mobile disloca- 
tions, on the other hand, give rise to plastic strain and 
dissipation. However they may not stay in their own 
category. When a mobile dislocation glides through the 
network of stored dislocations it can be pinned and be- 
come stored. This pinning can nevertheless be overcome 
by thermal activation and applied strain and the dislo- 
cation then becomes mobile again. 

A. Dislocation glide model 

In our dislocation glide model, we adopt a statistical 
representation in which individual dislocation character- 
istics (mobility, polarity, screw/edge character, line orien- 
tation, and position) are not treated explicitly. Of course, 
this is a great simplification of the complex problem of 
dislocations and presents the first step toward capturing 
the essential mechanisms and physics underlying dislo- 
cation glide in solid helium, similar to studies of other 
materials in the past. For simplicity, the large network 
of dislocations is represented as a continuous distribu- 
tion of dislocation lines within a representative volume, 
characterized by a dislocation density with units of total 
dislocation line length per volume. The total dislocation 
density in a material at time t consists of both mobile 
and stored (immobile) dislocations 

Ptot{t) = Pmob{t) + Pstored{t)- (1) 

In deformed crystals, stored dislocations assemble into 
networks within the interior, as shown in Fig.[TJ Here, we 
assume that stored dislocations take on a similar three 
dimensional pattern in strained solid helium. Accord- 
ingly, mobile dislocations moving in any direction and 
on any slip plane are likely to become immobile, when 
they encounter a forest of stored dislocations intersect- 
ing the glide plane, with which they strongly interact. 
When dislocations glide, it was proposed that the total 
dislocation density increases ]3 y35i36,38 

Aptotit) = Aep{t)/bL, (2) 

per unit time At. Here b is the magnitude of the Burgers 
vector, L is the mean- free path of mobile dislocations in 
a forest of dislocations, which is assumed to equal the 
interdislocation distance smd Aep is the incre- 

ment of plastic strain produced by dislocation glide. The 
increments in stored and mobile dislocation densities are 



generally proportional to Aptot- As portions of an ex- 
panding dislocation loop can react with other loops and 
become immobilized it changes the density of stored dis- 
locations by Ap stored- In this study, we assume for sim- 
plicity that newly added dislocations have equal chance 
of becoming mobile or stored. Thus, we arrive at the 
change per unit time At: Apmob = Apstored = Aptot/2. 

Mobile dislocations glide with a velocity distribution, 
which we parameterize only with its average speed v 
for convenience. In our dislocation network picture, 
see Fig. [TJ a dislocation can be pinned with potential 
Ui at local points, where it encounters other disloca- 
tions. We assume that the dislocation overcomes the 
pinning potential and glides by the combined effect of 
thermal activation (Arrhcnius term) and external drive 
Tapp- This assumption is supported by the experimen- 
tal observation that mechanical stress increases with ap- 
plied strain and decreases with temperaturci^i^'^ As 
the dislocation glides, the whole line still undergoes re- 
sistance from the Peierls stress Tpeieris and the mean-field 
dislocation-dislocation interaction stress Tint- The aver- 
age speed is thus described by a combination of thermo- 
dynamic potential and attempt frequency for overcoming 
the barrieri^2i^ 

v{t) = bi^DCxp[~AG{t)/kBT], (3) 

AG{t) ^ U^- [rappit) - Tpe.erls - Tint{t)] V, (4) 

where vu is the Debye or attempt frequency and V = b^A 
is the activation volume. A is assumed to be the average 
distance between pinning points, which equals the inter- 
dislocation distance L, and ks is the Boltzmann con- 
stant. The product bi'n is a typical sound speed of the 
solid, and AG is the thermodynamic Gibbs free energy 
for thermal activation of glide. Wc model the interac- 
tion between dislocations to be long ranged. In mean- 
field theory the interaction is inversely proportional to 
the interdislocation distance of a randomly oriented dis- 
location forest, which results in the resistive interaction 
stress22i32^d8 

Tint{t) = ap,by/ptot{t), (5) 

with elastic shear modulus p, and dislocation interaction 
coefficient a. Indeed, the dependence of Tint on ^ ptot 
has been observed in many materials i^i^^i^ 

We can now relate the macroscopic strain quantity 
with the above mentioned microscopic quantities of dis- 
location motion. The total applied strain is composed of 
both elastic and plastic strain, e^pp = fe + Cp, while the 
plastic strain rate due to glide obeys the kinematic equa- 
tion relating microscopic and macroscopic quantities^i^ 

tp(t) ^bv[t) p„ioh{t), (6) 

with the increment of plastic strain per time step At: 
Atp(t) = ep{t)At. The corresponding applied shear stress 
Tapp is given by Hookc's law of clastic deformation, 

Tapp{t) = pee{t) = p[eapp{t) - ep{t)] . (7) 
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This completes the set of thermodynamic and kinetic 
equations of our minimal dislocation glide model devel- 
oped for solid *He. 



B. Computation and model parameters 

The coupled system of equations are solved numeri- 
cally with a simple time-step evolution method, because 
it is sufficient to employ the Euler forward method to 
obtain the time evolution of the dynamics. In this work, 
we assume initial conditions eapp(O) = ep(0) = and 
Pmob(O) = Pstored{0) = 0.5ptot(0). First wc solve Eq. (7) 
for the increment of applied shear stress Arappit) at 
each time step At for an increment of applied strain 
Aeapp(t) ~ iappAt with given dc strain rate tapp, based 
on the difference between applied strain Capp and plastic 
strain €p. Then, we update the total applied shear stress 
Tapp{t) to get the average speed v{t) of dislocations from 
Eqs. (3) and (4). This way we solve the differential equa- 
tion for the plastic strain rate ep{t) in Eq. (6), which 
implicitly depends on e^, Tapp and Tint- Finally, we up- 
date the increment of the total dislocation density ptot (t) 
in Eq. (2) and the interaction stress Tint{t) in Eq. (5). 

In the calculations shown hereafter the pinning poten- 
tial Ui is set by the melting temperature = 1.86 K 
with the assumption that the pinning is from dislocation 
crossing or local melting. This strongly reduced value 
takes into account the quantum fluctuations, because in 
classical crystals Ui is of the order of the Debye temper- 
ature, which is roughly 26 K for hep ^He. For simplic- 
ity, we further assume that the Burgers vector is of the 
same magnitude as the lattice constant, b = 0.364 nm, 
which would correspond to glide on either the basal or 
prismatic planes of the hep crystal. We approximate the 
Debye frequency as i/d = Qoks/h = 600 GHz with a 
corresponding sound speed of bvo = 218 m/s and shear 
modulus p = 13.7 MPa.— For the dislocation interac- 
tion coefficient, we use a = 0.01. This value is much 
smaller than for the hep metal Zr, 0.1 < a < 2^^^ or 
fee Cu, 0.5 ^ a ^ 1 at room temperature , ^" i ■^^ A strongly 
reduced value should be expected when quantum fluctu- 
ations and interactions are important. In addition, it 
takes into account dislocation screening effects as well as 
the screening due to '^He atoms bound to the dislocation 
line. 

The Peierls barrier is inferred from the recent ac shear 
displacement experiments in which the softening hap- 
pens at an applied strain larger than 2.0 x 10~®fii^i^ 
By interpreting the softening of the solid as disloca- 
tion glide, we obtain Tp^ieris ~ 1.0 x 10~^/i. This 
value is consistent with the experiment by Sanders and 
coworkers^ who reported a nearly vanishing Peierls 
stress, whereas Hiki and Tsuruoka— reported a larger, 
but still small Tp^ieris ~ 10~^/x. On the other hand, 
Monte Carlo simulations predicted much larger values 
for Tpeieris ~ 0.006 — 0.03^ for edge and screw disloca- 
tions, respectively)^ Clearly the experimentally reported 



values for the Peierls stress are unusually small, which 
might be the combined effects of quantum fluctuations, 
which play a role at low temperatures in reducing the 
barrier for the escape of immobile dislocations and the 
nucleation of double kinks or small loops of atomic size 
6 vs. A 1^ If in fact the activation volume for overcoming 
the Peierls barrier is on the order of V = rather than 
V = 6^A, then the Peierls stress was effectively reduced 
by approximately 6/ A in our formulation of the energy 
barrier AG in Eq. (4). Everything else the same, we can 
either write for the potential energy of the Peierls barrier 
Ep = Tpeieris X 6^ A Or Ep = Tp{{^^i^ X 6^. Consequently, 
^^eris - Tpe^erls{A/b) = IQ-'^pA/b - 0.03^ for typ- 
ical dislocation density ptot = l/A'^ — 10^ m^^. This 

estimate of T^/ieris ^1^° i'^ reasonable agreement with 
the Monte Carlo simulationsi^ Thus a systematic study 
of the Peierls stress or critical strain versus dislocation 
density may resolve the true value of the Peierls stress. 

Of course, all material parameters used above arc sub- 
ject to updates in experiments. We believe that the most 
realistic result can be obtained once the material param- 
eters can be verified. In the next section, we will study 
how the plastic nature is revealed in the WHR and other 
experimental observables as well as a systematic study of 
the parameter dependence on the results. 



III. RESULTS AND DISCUSSION 
A. Strain-stress curves 

Dislocation glide manifests itself macroscopically as a 
change in the WHR. Prior to the onset of dislocation 
glide, the solid behaves elastically and the WHR equals 
the elastic shear modulus. At sufficiently large applied 
deformation, some of the pre-existing stored dislocations 
become unpinned and glide in reaction to the external 
stress. The glide of dislocations accommodates the defor- 
mation strain and the slope of the stress-strain response 
deviates from the shear modulus noticeably. The point 
of departure on the stress-strain response is normally re- 
ferred to as the yield point, with corresponding yield 
strain and yield stress. Above the yield point is the so- 
called plastic regime. In this study, the yield point of the 
stress-strain curve is de&ied as the point where the ratio 
of applied stress over applied strain equals 0.99//. Al- 
though this criterion is arbitrary, it is easy to implement 
and does not change any of the key results or conclusions. 

The calculated results are in qualitative agreement 
with recent ac shear experimentSfi'^'^ In Fig.[2Ua) we plot 
the predicted shear stress-strain curves for solid ^He de- 
formed under different temperatures with an initial dis- 
location density po = 1-0 x 10^ m~^i^ We ffiid that yield 
strains at all temperatures are roughly the same value of 
2.6 X 10^^, which is close to the critical strain 2.2 x 10~^ 
in the ac experiments. 

The temperature effect is most apparent in the stress- 
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FIG. 2: (color online), (a) Stress-strain curves for different temperatures at an applied strain rate of eapp = 10"'' s~^ with 
initial total and mobile dislocation densities of ptot = 1.0 x 10^ va"^ and pmob = 0.5 x l(f m~'^ . (b) Plot of the corresponding 
average speed-strain curves of (a) are nearly identical at different temperatures, (c) Plot of the corresponding dislocation 
density-strain curves of (a) are nearly identical at different temperatures, (d) The "effective" shear modulus /ie// = Tapp/eapp 
versus temperature; data for different strains have been shifted to coincide at 300 mK. 



strain curve above the yield stress. At high temperature 
the WHR approaches zero, which means that the mate- 
rial undergoes great plastic deformation even with mini- 
mum applied stress. In Fig.[51Jb) and (c), the temperature 
effect on the evolution of the average speed and disloca- 
tion density are substantially weaker than the ones shown 
for the parameters in panel (a) . There are two reasons for 
this difference: (1) The dislocation density in panel (c) 
is at the level of 10^ m^^; thus, a small difference on the 
order of 10^ m~^ still has a strong influence on the plastic 
strain rate calculated from Eq. (6), especially when the 
applied strain rate is iapp = 10~^ s~^. (2) The difference 
in stress-strain curves in Fig.[2l^a) is a consequence of the 
accumulated difference of plastic strain rates from each 
time step. Thus, the temperature effect is much more 
prominent in the stress-strain curves. 

One can visualize the microscopic processes of the de- 
formation through macroscopic quantities. The average 
dislocation speed v is of the order of ~ 2 m/s or 1% 
of the speed of sound and is plotted in Fig. [DJb). It 
shows that below the yield point v is zero and a boost 
is found right at the yield point. A slow decrease in v 
above the yield strain occurs because the increasing dis- 
location density increases the mean-field interaction and 
sets a higher effective barrier for gliding. Notice that al- 



though the "average" dislocation speed decreases, there 
are still fast dislocations to accommodate the increasing 
external strain; this has been seen in an earlier study us- 
ing discrete dislocation dynamicsi^ On the other hand, 
the increase in dislocation density p, see Fig. [2][c), is not 
as dramatic as for w, which we found to be true for all of 
our calculations. 

As seen in Fig. [2ja), the model also predicts a strong 
temperature dependence for the stress-strain curve, 
which follows directly from the kinetics embedded in the 
model of Eq. (2). The rate of dislocation motion across 
the solid is assumed to be controlled by thermally acti- 
vated glide, Eq. (2), and therefore, higher temperatures 
facilitate dislocation glide, enabling dislocations to glide 
faster at lower stresses. 

To summarize and allow direct comparison with exper- 
iments, we re-plot Fig. HJa) as an effective shear mod- 
ulus /ie// = Tapp/^app in Fig. [2jd). Note that this is 
the definition used in the ac shear experiments, rather 
than the WHR, irrespective of the value of the applied 
strain. One can see that the /^e// curve decreases with 
increasing temperature, while the zero-temperature value 
/ie//(0) decreases with increasing applied strain, which is 
in qualitative accordance with experiments. It is encour- 
aging that our dc strain-rate model actually captures the 
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essential features of the ac shear modulus experiments. 

B. Dislocation density dependence 

We show through varying the initial dislocation 
density that our model describes reported annealing 
experiments.— Dislocations are created during the growth 
process and can either be increased or decreased by pre- 
stressing or annealing. In the calculation we focus on 
two initial dislocation line densities of high density (HD), 
2.0x10*^ m-^ and low density (LD), 1.0x10^ m'^, which 
correspond to an as-grown and annealed sample, respec- 
tively. We see again the plastic behavior above the yield 
point. The yield point is higher for the HD case, see 
Fig.ElJa). This can be understood since the larger mean- 
field interaction in the HD sample possesses a stronger 
resistance to gliding, because with increasing density the 
interaction stress Tint increases as described by Eq. (3). 
Notice that the increase in dislocation density after ap- 
plying external strain is negligible, as shown in Fig.[3jb). 
This behavior is thus mainly determined by the initial 
dislocation density. In general the HD case also shows a 
larger crossover zone from the elastic to plastic regime, 
since a larger activation is required for dislocation lines 
to glide. We also compare the temperature effect on the 
two different density samples. In the temperature range 
of interest, 20-200 mK, the distinction between the two 
densities is clear. To better see the temperature effect, 
we plot the WHR as a function of temperature (Fig.[3l[d)) 
at a given applied strain, Capp = 2.65 x 10^^. The WHR 
of the HD case deviates from the zero-temperature shear 
modulus at higher temperatures, approximately at 100 
mK, compared to the LD scenario, which deviates at 
approximately 20 mK. This is because the larger dislo- 
cation interaction sets a higher barrier that requires a 
higher temperature to thermally activate glide. Again 
it is encouraging to see that our simple dislocation glide 
model captures the essential features of the ac shear ex- 
periments between as-grown and annealed samplesi^ 

C. Dissipative glide 

Dissipation is one of the essential features of disloca- 
tion glide described in this paper. To overcome barriers 
and obstacles when gliding, a dislocation dissipates en- 
ergy Wp = Tapp Aep per time step At. The corresponding 
dissipation rate is expressed as Wp . From this expression 
it is easy to see that significant dissipation or damping 
happens only in the plastic regime above the yield point 
Ty, where Cp > 0. In general, a dissipation peak accompa- 
nies a drastic change in the WHR as a function of temper- 
ature. The plastic dissipation corresponding to the WHR 
in Fig. [3jd) is shown in Fig. |4l Here, the dissipation peak 
is centered near the temperature, where the WHR, dr/de, 
is changing most rapidly, see Fig. [31^d). This behavior is 
similar to the one seen in the inverse quality factor 1/Q 



of ac shear experiments. In addition, the position of the 
peak shifts to lower temperatures as the dislocation den- 
sity decreases from the HD to LD scenario and the height 
of the peak increases with decreasing dislocation density. 
This narrowing of the dissipation peak can be understood 
in terms of fewer dislocations in the system resulting in 
fewer mobile dislocations with larger mean-free path to 
move. Again, that in turn causes a faster release of stored 
elastic energy in the system. 



D. Strain rate dependence 

Next, we study the effect of the key control param- 
eter in experiment on the plastic response of solid ^He, 
namely, the applied strain rate iapp ■ Thermally activated 
rate-controlled dislocation glide implies a sensitivity to 
the applied strain rate. Two different strain rates were 
used in the following calculations: a low strain rate (LSR) 
of 10"'' s-i and a high strain rate (HSR) of lO'^ s-\ In 
Fig. [5lja)-(c) we compare the model results for low and 
high strain rates. From Fig.[Sl[a) we can see that increas- 
ing the strain rate increases the yield stress and yield 
strain. In other words, as the strain rate increases, the 
transition from the elastic to plastic regime is delayed. 

This delay is due to the failure of plastic deformation to 
respond fast enough (at each time step) to the large ap- 
plied strain rate for HSR. The plastic strain rate, ip oc v, 
is roughly ten times larger for HSR (Fig. [5] (b)) while 
the applied strain is in fact a hundred times larger. The 
plastic deformation, therefore, does not affect the ap- 
plied stress much until a much larger strain is applied, 
see Eq. ([7]); the yield point then happens at a higher 
drive. 

Further, we observe that the sensitivity of the yield 
stress and strain to strain rate is enhanced at higher tem- 
peratures. Since the yield stress in the high strain rate 
case is higher, the driving force for dislocation glide is 
larger, causing the dislocations on average to speed up 
as shown in Fig. ^h) . This sensitivity can also be seen 
in the variation in the WHR, dr/de, shown in Fig. [SKc). 
Here we see that the extent of the elastic regime is wider 
at the higher temperature. Fig.[Sl[d) compares the dr/de 
from high to low temperatures at different strain rates. 
Two features are revealed. First, changing the strain rate 
does not affect the low temperature behavior, i.e., below 
20 mK. Second, compared to the HSR case, the WHR 
decreases more rapidly in the LSR case, resulting in a 
larger drop of the WHR at the same amount of applied 
strain. Notably, the latter feature can also be found in 
the ac shear experimentsjii^ although the impact on the 
WHR is not as obvious as in our calculations. The differ- 
ence may be related to the usage of fi^f f rather than the 
WHR in experiments or the different loading methods in 
ac versus dc experiments or the boundary conditions4^ 
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FIG. 3: (color online), (a) Stress-strain curves for high (HD) and low (LD) initial dislocation densities at ta 



10" 



(b) 



Plot of the corresponding dislocation density-strain curves of (a), (c) Plot of the corresponding WHR-strain curves of (a), (d) 
The WHR versus T for HD and LD cases. The curves for HD and LD cases correspond to the initial total dislocation density 



of 1.0 X 10" m"^ and 2.0 x 10" m"^, respectively. AH data are taken from the same strain amplitude at ea 



2.650 X 10" 
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FIG. 4: (color online). Plastic dissipation curves correspond- 
ing to the WHR curves of Fig. [Sjjd) for low (LD) and high 
(HD) densities of dislocations. 



E. Model parameter dependence 

Finally, it is important to characterize the material 
parameter dependence of our model results. A suitable 



experimentally observable quantity for that purpose is 
the yield point. We study systematically the effect of the 
pinning potential Ui, the Peierls barrier Tpeieris, and the 
dislocation interaction strength a at 20 mK and 300 mK, 
since they play important roles in determining the aver- 
age dislocation speed v given in Eq. (3) and the WHR. 
We keep all model parameters the same as in Sec. III.A, 
except for the parameter whose dependence on the yield 
stress is investigated. 

All parameters represent physical barriers to disloca- 
tion glide. Accordingly, the yield stress Ty increases with 
all three parameters as expected (Fig. IHl El ©■ Among 
the three parameters, the dependence on Ui is least cru- 
cial, see Fig. [SJ since the value of Ui is known to be 
restricted by the melting temperature Tm = 1-86 K. (1) 
The yield stress Ty varies by less than an order of mag- 
nitude in the given range of Ui values with the scaling 
relation Ty ^ O.lUi^. (2) The interaction strength 
a, although difficult to control in experiments, is be- 
lieved to range between 0.1 to 2 when compared with 
reported values for ordinary metals. In our model calcu- 
lations we chose a = 0.01 to properly describe the WHR. 
This strongly reduced value, compared to ordinary met- 
als, might be caused by interdislocation screening effects 
originating from quantum fluctuations. Our calculation 
shows that Ty varies by roughly an order of magnitude 
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FIG. 5: (color online), (a) Stress-strain curves for different applied strain rates at iapp = 10^^ s^^; (b) plot of the corresponding 
average speed-strain curves of (a); (c) plot of the corresponding work-hardening rate-strain curves of (a); (d) work-hardening 
rate, WHR, versus temperature for different applied strain rates. (The curves for high and low applied strain rates correspond 
to 10"'^ and 10"'* s~^, respectively). All data are taken from the same strain amplitude at eapp ~ 2.650 x 10~*. 
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FIG. 6: (color online). Plot of yield stress Ty versus pinning 
potential Ui at both low (20 mK) and high (300 mK) temper- 
atures. 



FIG. 7: (color online). Plot of yield stress Ty versus dislo- 
cation interaction strength a at both low (20 mK) and high 
(300 mK) temperatures. 



for a between 0.01 and 1, see Fig. [T] The curve satu- 
rates befow a < 0.01, because then Ty is dominated by 
other barrier types, e.g., Peierls barrier. For a > 10~^ 
an approximate scaUng of Ty ~ (0.14 -I- 5a) Pa can be 



inferred. (3) For the Peierls potential, however, the yield 
stress changes by more than three orders of magnitude 
over the reported range shown in Fig. [51 It satisfies a sim- 
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FIG. 8: (color online). Plot of yield stress Ty versus pinning 
potential rpeieria at both low (20 mK) and high (300 mK) 
temperatures. 

pie scaling relation OC 1.4 • lO'^ Tpeierls/ fJ-Pa ~ Tpeierls, 

which can be used to estimate the order of magnitude 
of the Peierls barrier. From this work it has become 
clear that a reliable and independent determination of 
the Peierls barrier is needed. On the other hand, increas- 
ing temperature accentuates the parameter dependence 
on Ty in general, because the solid is softer at higher tem- 
peratures. Fortunately, the temperature dependence of 
Ty versus Tpeieris , the least well-known parameter in our 
dislocation model, is negligible. 



solid helium. No supersolidity is required to explain these 
effects. 

One may speculate about the realm of quantum and 
super plasticity in the quantum solid helium. So far 
the results of our thermally activated dislocation glide 
model would suggest that quantum tunneling of disloca- 
tion lines has not yet been seen in experiments. Indeed, 
measurements of the dynamic shear modulus down to 20 
mK confirm the dominance of thermally activated relax- 
ation behavior in the dissipation processj^i^i^ 

Our predictions for the dc strain rate measurements are 
obtained with very reasonable parameter values. Since 
the parameter range is generally constrained, the model 
gives rather robust predictions, which is encouraging par- 
ticularly in the event that some of the parameters need 
to be updated with the advent of more accurate materi- 
als characterization. Future dc experiments can certainly 
narrow down the range of values for some of the model 
parameters. For example, our model calculations show 
a simple linear relation between yield stress and Peierls 

stress, Ty ~ ^Peierls- 

Although the current dislocation glide model does not 
predict hysteresis, it can be extended, in principle, to 
capture the irreversible processes of creation and annihi- 
lation of dislocations during loading and reverse loading. 
We encourage and welcome experiments that will allow a 
more precise dynamical and structural characterization 
of the quantum solid '^He to provide stringent tests on 
quantum vs. classical dislocation glide at low tempera- 
tures. 



IV. CONCLUSION 

In summary, we proposed and studied a model of dc 
strain rate driven dislocation glide in solid ^He. We pre- 
dict stress-strain curves and work-hardening rates as well 
as the role of the dislocation interactions and Peierls 
barriers. We find that the shear modulus can exhibit, 
what appears to be, anomalous elasticity. However, these 
anomalies can be described rather well by classical theo- 
ries of plastic deformation with renormalized model pa- 
rameters due to zero-point fluctuations of the quantum 
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samples. Application was made to single-crystal fee Ni 
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